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74 Abstract 

75 The OPERA neutrino detector in the underground Gran Sasso Laboratory (LNGS) was designed 

76 to perform the first detection of neutrino oscillations in appearance mode through the study of z/^ Vt 

77 oscillations. The apparatus consists of an emulsion/lead target complemented by electronic detectors 

78 and it is placed in the high energy long-baseline CERN to LNGS beam (CNGS) 730 km away from 

79 the neutrino source. Runs with CNGS neutrinos were successfully carried out in 2007 and 2008 with 

80 the detector fully operational with its related facilities for the emulsion handling and analysis. After a 

81 brief description of the beam and of the experimental setup we report on the collection, reconstruction 

82 and analysis procedures of first samples of neutrino interaction events. 



83 1 Introduction 

84 Neutrino oscillations were anticipated nearly 50 years ago [T] but they have been unambiguously observed only 

85 recently. Several experiments carried out in the last decades with atmospheric and accelerator neutrinos, as well 

86 as with solar and reactor neutrinos, contributed to our present understanding of neutrino mixing (see e.g. [2] for a 

87 review) . 
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Figure 1: View of the OPERA detector. The upper horizontal hnes indicate the position of the two 
identical supermodules (SMI and SM2). The "target area" is made of walls filled with ECC bricks 
interleaved with planes of plastic scintillators (TT). Arrows show the position of the VETO planes, the 
drift tubes (PT) pulled alongside the XPC, the magnets and the RPC installed between the magnet iron 
slabs. The Brick Manipulator System (BMS) is also visible. See [6] for more details. 



B8 As far as the atmospheric neutrino sector is concerned, accelerator experiments can probe the same oscillation 

89 parameter region as atmospheric neutrino experiments [3] . This is the case of the OPERA experiment that has the 

90 main scientific task of the first direct detection of j/^ Vt appearance, an important missing tile in the oscillation 

91 scenario [4l[5l[6]. 

92 OPERA uses the long-baseUne (L=730 km) CNGS neutrino beam [7] from CERN to LNGS, the largest un- 

93 derground physics laboratory in the world. The challenge of the experiment is to measure the appearance of 

94 from i/^ oscillations in an almost pure muon-neutrino beam. Therefore, the detection of the short-lived r lepton 

95 (cr — 87.11 /im) produced in the charged-current (CC) interaction of a is mandatory. This sets two conflicting 

96 requirements: a large target mass to collect enough statistics and an extremely high spatial accuracy to observe 

97 the short-lived t lepton. 

98 The T is identified by the detection of its characteristic decay topologies either in one prong (electron, muon 

99 or hadron) or in three- prongs; its short track is measured with a large mass target made of 1 mm thick lead plates 

100 (target mass and absorber material) interspaced with thin nuclear emulsion films (high- accuracy tracking devices). 

101 This detector is historically called Emulsion Cloud Chamber (ECC). Among past applications it was successfully 

102 used in the DONUT experiment for the first direct observation of the [Si . 

103 OPERA is a hybrid detector made of two identical Super Modules (SM) each consisting of a target section of 

104 about 625 tons made of emulsion/lead ECC modules (hereafter called "bricks"), of a scintillator tracker detector 

105 (TT) needed to trigger the read-out and localize neutrino interactions within the target, and of a muon spectrometer 

106 (Figure[l]). The detector is equipped with an automatic machine (the Brick Manipulator System, BMS) that allows 
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Figure 2: Schematic view of two bricks with their Changeable Sheets and target tracker planes. 



107 the online removal of bricks from the detector. Ancillary, large facilities are used for the handling, the development 

108 and the scanning of the emulsion films. Emulsion scanning is performed with two different types of automatic 

109 microscopes: the European Scanning System (ESS) [21 [TO] and the Japanese S-UTS [TT]. 

no A target brick consists of 56 lead plates of 1 mm thickness interleaved with 57 emulsion films [T^]- The plate 

111 material is a lead alloy with a small calcium content to improve its mechanical properties [13j . The transverse 

112 dimensions of a brick are 12.8 x 10.2 cm^ and the thickness along the beam direction is 7.9 cm (about 10 radiation 

113 lengths). The bricks are housed in a light support structure placed between consecutive TT walls. More details on 

114 the detector and on the ancillary facilities are given in [B] . 

115 In order to reduce the emulsion scanning load the use of Changeable Sheets (CS) film interfaces [2] , successfully 

116 applied in the CHORUS experiment 15J, was extended to OPERA. Tightly packed doublets of emulsion films are 

117 glued to the downstream face of each brick and can be removed without opening the brick. The global layout of 

118 brick, CS and TT is schematically shown in Figure [H 

119 Charged particles from a neutrino interaction in the brick cross the CS and produce a signal in the TT 

120 scintillators. The corresponding brick is then extracted and the CS developed and analyzed in the scanning 

121 facilities at LNGS and in Nagoya. The information of the CS is then used for a precise prediction of the position 

122 of the tracks in the most downstream films of the brick, hence guiding the scan-back vertex finding procedure. 

123 A reconstructed CC event is shown in the bottom panels of Figure[3] In this case the detached event dimensions 

124 are of the order of a few millimeters, to be compared with the 10 m scale of the whole event reconstructed with 

125 the electronic detectors (top panels of Figure [3]) . 

126 First neutrino data were collected by OPERA in 2006 [5] with the electronic detectors alone, and then in 

127 2007 and 2008, for the first time with target bricks installed. All steps from the prediction of the brick where the 

128 interaction occurred down to the kinematical analysis of the neutrino interactions are described in the following 

129 using as benchmark a sub-sample of the statistics accumulated during the CNGS runs. The procedure has proven to 

130 be successful. We are presently in the process of a quantitative evaluation of the different experimental efficiencies 

131 that are involved in the analysis procedure, profiting of the increasing statistics of the reconstructed neutrino 

132 events. 
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Figure 3: Top panels: on line display of an event seen by the OPERA electronic detectors (side and top 
views): a interacts in one of the first bricks of the first supermodule (SM) yielding hadrons and a 
muon which is detected in both SMs and whose momentum is measured by the magnets of the two SMs. 
Bottom panels: the vertex of the same event observed in the emulsion films (side, top and front views). 
Note the two 7 e^e~ vertices: the opening angle between them is about 300 mrad. By measuring the 
energy of the 7's one obtains a reconstructed invariant mass of 110 it 30 MeV/cP, consistent with the tt" 
mass. 



5 




Date 



Figure 4: Integrated number of protons on target (p.o.t.) as a function of time for the 2008 CNGS run 
(June-November) . 



133 2 Real time detection of the CNGS Beam 

134 The CNGS neutrino beam [7| was designed and optimized for the study of i/^ — > oscihations in appearance 

135 mode by maximizing the number of CC v-r interactions at the LNGS site. After a short commissioning run in 

136 2006 the CNGS operation started on September 2007 at rather low intensity. The first event inside the OPERA 

137 target was observed on October 3'"'^. Unfortunately, due to a fault of the CNGS facility, the physics run lasted 

138 only a few days. During this run 0.082 x 10^^ protons on target (p.o.t.) were accumulated with a mean value of 

139 1.8 X 10^'^ protons per extractiorQ: this corresponds to about ^ 3.6 effective nominal days of running. With such 

140 an integrated intensity 32 neutrino interactions in the bricks and 3 in the scintillator material of the target tracker 

141 were expected; we actually observed 38 events on time with the arrival of the beam at Gran Sasso. 

142 A much longer run took place in 2008 when 1.782 x 10^^ protons were delivered on the CNGS target. OPERA 

143 collected 10100 events on time and among them 1700 interactions in the bricks. The other events originated outside 

144 the target region (spectrometers, supporting structures, rock surrounding the cavern, hall structures, etc.). The 

145 run featured a poor initial efficiency of the CERN complex, with an average of about 40% until reaching an average 

146 value of about 60%. From October to November OPERA gathered the same number of events as from June to 

147 August. The 2008 CNGS integrated p.o.t. intensity as a function of time is shown in Figured 

148 During the 2007 and 2008 runs all electronic detectors were operational and the live time of the data acquisition 

149 system exceeded 99%. More than 10 million events were collected by applying a minimum bias filter. The selection 

150 of beam related events relies upon a time stamp, based on the time synchronization accuracy of 100 ns between 

151 the CERN beam GPS tagging and the OPERA timing system. 

152 An automatic classification algorithm provides high efficiency in the selection of neutrino interactions inside the 

153 OPERA target both for CC and neutral-current (NC) events at the expenses of a slight contamination of neutrino 

154 interactions in the external material. 

155 For the early 2007 run the algorithm selected 53 events occurring inside the target while expecting 50; the 

156 contamination from neutrino interactions outside the target amounted to 37% (Monte Carlo simulation). The low 

157 purity of the selected event sample was due to some sub-detectors still being in the commissioning phase and to the 

158 OPERA target that was only partially filled with bricks. In the 2008 run 1663 events were classified as interactions 

159 in the target (expected 1723) with a contamination from outside events of only 7%. 

160 The muon momentum distribution for events classified as CC interactions in the target is shown in the left 

161 panel of Figure [51 The distribution of the muon angle with respect to the horizontal axis is shown in the right 

^The 400 GeV proton beam is extracted from the CERN SPS in two 10.5 /is pulses, with design intensity of 2.4 x 10^^ 
p.o.t. 
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Figure 5: Left: momentum distribution of muons produced in CC neutrino interactions inside the OPERA 
target. Right: angular distribution of the muon tracks with respect to the horizontal axis. 



162 panel of Figure [SJ the beam direction angle is found to be tilted by 58 mrad, as expected from geodesy. 

163 An extensive study of the beam monitoring is being performed by using neutrino interactions both in the whole 

164 OPERA detector and in the surrounding rock material. This will be the subject of a forthcoming publication. 

165 3 Combined analysis of electronic detectors and nuclear emulsion film 

166 data 

167 We describe in the following the breakdown of the different steps carried out to analyze neutrino interaction events 

168 from the identification of the "fired" brick up to the detailed kinematical analysis of the vertex in the emulsion 

169 films. 

170 Once a trigger in the electronic detectors is selected to be compatible with an interaction inside a brick the 

171 following procedure is applied [6]: 

172 1. electronic detector data are processed by a software reconstruction program that selects the brick with the 

173 highest probability to contain the neutrino interaction vertex; 

174 2. this brick is removed from the target wall by the BMS and exposed to X-rays for film-to-film alignment. 

175 There are two independent X-ray exposures: the first one ensures a common reference system to the CS film 

176 doublet and the most downstream film of the brick (frontal exposure); the second one produces thick lateral 

177 marks on the brick edges, used for internal alignment and film numbering within the brick; 

178 3. after the first X-ray exposure the CS doublet is detached from the brick and developed underground, while 

179 the brick is kept in a box made of 5 cm thick iron shielding to reduce the radioactivity background; 

180 4. if the CS scanning detects tracks compatible with those reconstructed in the electronic detectors the second 

181 X-ray exposure (lateral marking) is performed and the brick is brought to the surface laboratory. The brick 

182 is then exposed to cosmic-rays for about 24 hours in a dedicated pit in order to select high-energy cosmic 

183 muons to provide straight tracks for a refined (sub-micrometric) film-to-film alignment; 

184 5. the brick emulsion films are then developed and dispatched to the various scanning laboratories in Europe 

185 and Japan. 

186 The procedure described above has proven to be successful. As an example, in Figure [5] we show the number of 

187 bricks extracted by the BMS per week, about one hundred. This is matched by the 100 CS developed and scanned 

188 per week in the LNGS (Italy) and Tono (Japan) scanning stations. 
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Figure 6: The number of bricks extracted per week by the BMS in 2008. 



1B9 Brick Finding and Changeable Sheet Interplay 

190 The efficiency for selecting the "fired" brick is the convolution of several effects and measurements. Here we 

191 discuss the two most important ones, the Brick Finding procedure and the Changeable Sheet measurement, for 

192 which preliminary results have been obtained from the analysis of partial samples of already scanned events. 

193 The brick finding algorithm exploits the tracking capabilities of the OPERA electronic detectors and, by 

194 combining this information with the output of a Neural Network for the selection of the most probable wall where 

195 the interaction occurred, provides a list of bricks with the associated probability that the interaction occurred 

196 therein. A preliminary estimate of the brick finding efficiency, limited to the extraction of the first most probable 

197 brick (for about 700 events) and not considering the small fraction (< 5%) of events for which the present electronic 

198 detector reconstruction fails, is compatible with the Monte Carlo estimate of 70% computed for a standard mixture 

199 of CC and NC events. A higher efficiency can be obtained by extracting also bricks ranked with lower probabilities. 

200 The tracking efficiency of single emulsion films can be measured by an exposure to high-energy pion beams 

201 and amounts to about 90% |10j. However, the measurement of the CS doublet efficiency in situ, in the OPERA 

202 detector, is by far more challenging, given the coarse resolution in the extrapolation of tracks from the electronic 

203 detectors to the CS. 

204 At present, we are studying the CS tracking efficiency by two independent approaches: (a) all tracks produced 

205 in already located neutrino vertices are followed downstream and searched for in the corresponding CS doublet; 

206 (b) muon tracks reconstructed by the electronic detectors and found in the CS are properly normalized to the total 

207 number of CC events where at least one track (not necessarily the muon) is found in the CS. The two methods 

208 yield a preliminary efficiency for finding a track in both films of the CS doublet which is compatible with the 

209 conservative expectation of 90% on a single film TU]. The experimental efficiency has been evaluated on a sample 

210 of 100 events scanned in both the European and the Japanese laboratories. We are presently working in order to 

211 further increase this efficiency by employing more advanced analysis techniques. 

212 Vertex analysis 

213 All tracks measured in the CS are sought in the most downstream films of the brick and followed back until 

214 they are not found in three consecutive films. The stopping point is considered as the signature either for a primary 

215 or a secondary vertex. The vertex is then confirmed by scanning a volume with a transverse size of 1 cm^ for 11 

216 films in total, upstream and downstream of the stopping point (see Figure [7]). Preliminary estimates of the vertex 

217 location efficiency are in agreement with the Monte Carlo expectations of 90% and 80% for CC and NC events, 

218 respectively. This evaluation has been performed on a sample of 500 located events. 
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Figure 7: Schematic view of the volume scan performed around the stopping point of the track. 




Impact Parameter (|xm) Multiplicity 

Figure 8: Left panel: impact parameter distribution of the muon track in CC events with respect to the 
reconstructed vertices. Right panel: charged track multiplicity distribution of the events. 



219 The track impact parameter distribution of the muon in CC events with respect to the reconstructed vertex 

220 position and the event track multiplicity distribution are shown in Figure [8l As expected, the impact parameter 

221 distribution is peaked at zero and has a mean value of 2.5 /im. The multiplicity distribution shows the anticipated 

222 enhancements for even track numbers due to the preferred interaction of neutrinos with neutrons. 

223 As an example, in Figures [51 and [TUl we show a NC and a CC event, respectively, fully reconstructed in the 

224 brick. A very " peculiar" event is shown in Figure 1111 the neutrino interaction occurred in the bottom layer of an 

225 emulsion film. Therefore, the associated nuclear fragments (large angle heavy ionizing tracks) are also visible in 

226 the film containing the vertex. 

227 Decay topologies 

228 Charm production and decay topology events have a great importance in OPERA for two main reasons. On 

229 the one hand in order to certify the observation of r events one should prove the ability of observing charm events 

230 at the correct expected rate. On the other hand, since charm decays exhibit the same topology as r decays, they 

231 are a potential source of background if the muon at the primary vertex is not identified (see Figure [12]). Therefore, 

232 searching for charm-decays in events with the primary muon correctly identified provides a direct measurement of 

233 this background. 
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Figure 9: Top panels: online display of one NC event seen by the OPERA electronic detectors. The 
regions filled with bricks are highlighted. Bottom panels: the emulsion reconstruction is shown in the 
bottom panels: top view (bottom left), side view (bottom center), front view (bottom right). 
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Figure 10: Top panels: online display of one CC event seen by the OPERA electronic detectors. Bottom 
panels: the emulsion reconstruction is shown in the bottom panels: top view (bottom left), side view 
(bottom center), front view (bottom right). 
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Figure 11: Display of the OPERA electronic detector of a CC interaction, top and side views. The 
emulsion reconstruction is shown in the middle panels: top view (bottom left), side view (bottom center), 
frontal view (bottom right). Bottom left panel : picture of the interaction vertex as seen by the microscope 
CMOS camera. The nuclear fragments produced in the interaction are visible. Bottom right panel 
: picture taken about 200 micron far from the interaction vertex. The minimum ionizing particles 
produced in the interaction are indicated by a circle. The muon track is also indicated. 
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Figure 12: Schematic view of the charm and tau decay topologies. 



234 Charm decay topologies were searched for in the sample of located neutrino interactions. Two events with 

235 charm-like topologies were found. By using the neutrino-induced charm-production cross-section measured by the 

236 CHORUS experiment ^17) about 3 charged-charm decays are expected to be observed in the analyzed sample. 

237 The event in Figure [T3] has high track multiplicity at the primary vertex and one of the scan-back tracks shows 

238 a kink topology. The measured decay angle is 204 mrad and the flight length of the decaying particle is 3247 /xm. 

239 The decay occurred in the third lead plate downstream of the interaction plate. No large angle tracks are produced 

240 at the decay vertex. This allows to further rule out the hadronic interaction hypothesis. The muon track and the 

241 charm candidate track lie in a back-to-back configuration (At/) ~ 165°) as one would expect for charm production. 

242 The daughter momentum, measured by using the Multiple Coulomb Scattering technique is 3.9^J g GeV/c at the 

243 90% C.L. Therefore, at the 90% C.L. the transverse momentum ranges between 600 MeV/c and 1150 MeV/c, well 

244 above the cut of 250 MeV/c applied to reject hadronic decays. According to the FLUKA Monte Carlo [18] the 

245 probability that a hadron interaction mimics a charm-decay with transverse momentum larger than 600 MeV/c is 

246 only 4 X 10^"^. 

247 The second charm-like topology is shown in Figure 1141 A 4-prong primary vertex is observed originating at a 

248 depth of about 30 /im in the upstream lead plate. The charmed hadron track points to a 3-prong decay vertex 

249 located at a distance of 1150 /im from the primary vertex (200 /xm inside the lead). All tracks have a clear CS tag. 

250 The interaction occurs downstream in the brick and the tracks only cross four emulsion films and the CS doublet 

251 (the two most downstream hits in the figure). The muon track and the charm candidate track lie in a back-to-back 

252 configuration (A(/i ~ 150°). The relativistic 7 of the charmed parent has been roughly estimated as the inverse of 

253 the average angle in space that the daughter tracks form with it. This leads to a 7 value of about 8.6, implying a 

254 parent high momentum of ~ 16 GeV/c (assuming the mass). The momenta of the daughter tracks have also 

255 been measured by extracting the downstream brick and using the Multiple Coulomb Scattering technique. The 

256 measured values are pi=2.4to.6, P2=1.3to 3 and p3= 1.2lJ;4 GeV/c (transverse momenta of about 610, 90 and 

257 340 MeV/c, total momentum: 4.8j^Q'g GeV/c), at the 90% C.L. The probability of a hadron interaction has been 

258 evaluated using FLUKA and amounts to 10~^. Assuming a. D ^ Ktttt decay, an invariant mass of I.iI^q'^ GeV/c^ 

259 is obtained. On the other hand assuming a, Dg ^ KKir decay an invariant mass of 1.5^g^ GeV/c^ is derived. In 

260 the latter case the invariant mass is consistent with the mass of a charmed hadron while in the second case the 

261 consistency is marginal. The probability of a decay in flight oi a, K is about 10^^. 
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Figure 13: Online display of the OPERA electronic detector of a charged-current interaction with a 
charm- like topology (top panel). The emulsion reconstruction is shown in the bottom panels where the 
charm-like topology is seen as a track with a kink: top view (bottom left), side view (bottom center), 
frontal view (bottom right). 
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Figure 14: Online display of the OPERA electronic detector of a charged-current interaction with a 
charm-like topology (top panel). The emulsion reconstruction is shown in the bottom panels where the 
charm- like topology is seen as a three-prongs secondary vertex: top view (bottom left), side view (bottom 
center), frontal view (bottom right). 
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262 4 Conclusions 



263 The 2007 and 2008 CNGS runs constitute an important milestone for the LNGS OPERA experiment searching for 

264 Vr oscillations. First samples of neutrino interaction events have been collected in the emulsion/lead target 

265 and allowed to check the complete analysis chain starting from the trigger down to the neutrino vertex location in 

266 the emulsions and to the topological and kinematical characterization of the event. 

267 In this paper we reported on the capability in performing an online identification, extraction and development 

268 of the bricks where the neutrino interaction occurred, vertex location and kinematical reconstruction. The overall 

269 performance of the experiment during the running phase and through the analysis chain can be summarized by 

270 stating that: 

271 • all electronic detectors performed excellently allowing the precise localization of the brick hit by the neutrino; 

272 • the electronic detector event reconstruction was tuned to the brick finding procedure which operated for the 

273 first time with real neutrino events providing good results; 

274 • all experimental activities from brick removal upon identification to the X and cosmic-ray exposures, brick 

275 disassembly and emulsion development, have been successfully accomplished. At present more than 100 

276 bricks per week can be routinely handled; 

277 • the scanning of the Changeable Sheets can be performed with the expected detection efficiencies; 

278 • vertex location was successfully attempted for both CC and NC events; 

279 • the topological and kinematical analyses of the vertices were successfully exploited and led to an unambiguous 

280 interpretation of neutrino interactions. In particular two events with a charm-like topology were found so 

281 far in the analyzed sample. This is fully consistent with expectations based on the known neutrino-induced 

282 charm production cross-section. 

283 The above considerations make us confident that the OPERA experiment is definitely in its production phase 

284 with the CNGS beam and that the scene has been set for the discovery of the t appearance. 
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